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1. Introduction 

A quantum system in n-dimensions is integrable if it admits n integrals of motion that 
are in involution. It is superintegrable if it possesses additional constants of the motion, 
i.e. more than there are degrees of freedom. In such cases, only subsets of n of the 
symmetry generators can commute amongst themselves. A model is called maximally 
superintegrable if it has 2n — 1 conserved quantities, including the Hamiltonian, that 
are algebraically independent. 

In view of the special properties and applications that superintegrable models 
enjoy, there is considerable interest in enlarging the set of documented systems in 
that category. While not doing justice to the many contributors to this field, let 
us mention in this connection that the recent identification of the TTW family [1] 
in particular, has attracted a lot of attention. It stimulated much effort aimed at 
demonstrating the superintegrability of (the higher-order) members of the family [2], 
obtaining the constants of the motion and determining the structure of the symmetry 
algebras pi H]. To that end, the recurrence approach [5] proved most general and 
powerful for determining the explicit integrals of the quantum systems. 
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We here wish to add to this body of knowledge by introducing another infinite 
family of maximally superintegrable systems in the plane. This new family shares many 
features with the TTW one; it has however a distinctive trait in that its Hamiltonians 
involve a refiection operator. 

Hamiltonians with refiection symmetry are customary in the context of the 
Calogero-Sutherland integrable models and their generalizations. The exchange 
operator formalism is indeed best suited to express the conserved quantities of these 
systems, to show that they are in involution and to find the ladder operators [HI [TJ [8]. 
Interestingly, it is precisely this formalism that has been used in [2] to address the 
superintegrability of the TTW system. It should also be recalled that the creation and 
annihilation operators of the parabose oscillator are realized with refiection operators 
P, [ini HI] . This background provides reasons to study problems with refiection operators 
in addition to the intrinsic interest that finding new superintegrable models entails. 

The increasing appreciation that certain systems require the use of refiection 
operators for their description [12], [131 E] is not unrelated to the fact that classical 
orthogonal polynomials that are eigenfunctions of first-order differential-difference 
operators with reflections have recently been discovered [151 IBl 113 HE] • The simplest 
among those are the little -1 Jacobi polynomials which turn out to be g — t- — 1 limits 
of the little q- Jacobi polynomials [15] . The identification and characterization of these 
new classes of polynomials is likely to enlarge the set of problems that are deemed 
exactly-solvable. 

We have in fact already found this to be the case in a recent article [19j where we 
have considered a realization of supersymmetric quantum mechanics, that is achieved 
by using as supercharges such differential-difference operators with reflections. We have 
obtained through this approach an exactly solvable supersymmetric extension of the 
Scarf I Hamiltonian and found that its wave functions are given in terms of little -1 
Jacobi polynomials. As it turns out, this system will be a key building block in the 
construction of our superintegrable family since it will essentially be the angular part 
of our 2-d Hamiltonians and play the role of the constant of motion responsible for 
separation of variables in polar coordinates. As we soon shall see, the Hamiltonians of 
the superintegrable models in our set, are formed by combining the radial part of the 
2-d harmonic oscillator with this new, exactly- solvable, supersymmetric Hamiltonians 
in 1-d. We surmise that this illustrates a systematic way to construct problems that are 
exactly solvable and potentially superintegrable. 

The outline of the paper is as follows. In section 2, we introduce the family of 
Hamiltonians that is the object of this article. It comprises an inflnite number of systems 
labeled by a real constant k and further involves 3 real parameters a;, a and /3. The exact 
solutions of the associated Schrodinger equations are provided in section 3. As already 
indicated, the little -1 Jacobi polynomials will appear. For rational k, the energy spectra 
will show higher degeneracies, signaling the superintegrability of these systems which 
will be proven in section 4 using the recurrence approach and proven directly in the 
appendix. The (higher) integrals of motion will be explicitly obtained and the structure 
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relations of the symmetry algebra will be determined. In section 5, we shall indicate how 
another infinite family of superintegrable models can be obtained from the one under 
consideration by trading, via coupling constant metamorphosis, the harmonic potential 
for the Coulomb one in the radial part [20] . 



2. The infinite family of Hamiltonians 

We shall consider the following ensemble of Hamiltonians in M? written in polar 
coordinates (r, 9) : 

r 2 \ sm kO J 2 \r'^ cos"' kU I 

where R is the reflection operators with respect to 6 

Rf{0) = f{-e). (2.2) 

We require the following restrictions on the real parameters a, /3, and k: a > —1, /3 > — 1, 
k 0. We must also have r G [0, oo) and — ^ < 6* < ^. 

The Hamiltonians Hk{r,6]u,0, P) obtained by setting a = do not involve the 
reflection operator. The conclusions of our analysis will prevail also for this subfamily 
and, as in the case a 7^ 0, all these systems will be shown to be superintegrable. 

When k = 1, Hi has the following expression in Cartesian coordinates (xi = 
r cos 6', X2 = r sin 6') : 

= -S..-a.,+. + _ j , (2.3) 

where R2 is the reflection operator in X2- Now let 

TT TT 
- <(h< - 

2-^-2 



k9, --<</)<-, (2.4) 



we then have 



Hk = -dl - -dr + a;V + %H^, (2.5) 



where 

ak' / f-cos</>i? \ / f-sin0 \ 
- + — I, sin^^ J + — ) ■ ^^-^^ 

The angular operator corresponds (up to a factor) to the exactly solvable 1- 
d Hamiltonian that was studied in [19]. It is manifestly an extension of the original 
Scarf I Hamiltonian [21]. It was shown to be supersymmetric within a framework where 
supercharges are taken to be differential- difference operators with reflections [22[ [T9] . In 
the present case, it is readily checked that 

= Qlp (2.7) 

with 

= (d^ - R - (2.8) 

V 2 cos (hi 2 sm (B 
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Obviously, the Hamiltonians have been constructed by combining the radial part 
of the 2-d harmonic oscillator Schrodinger operator with this supersymmetric 1-d 
Hamiltonian. is thus set as a constant of motion and is therefore integrable (at 
a minimum). The exact solvability of the eigenvalue equations associated to ensures 
that the 2-d Hamiltonian also shares that property through separation of variables as 
we shall see next. 



3. Exact solvability 



We are looking for the solutions \E'(r, 0) to the Schrodinger equation 

i/fc^(r,0) = Evl/(r,0) (3.1) 

with Hk as in (12. 5p and (12.61) . (The solutions to the original problem are obtained by 
making the substitution (12.41) .) Let us separate the variables and write 



^(r,0) = fi(r)<l>(0) 

The separated equations read 

1 „ ka^ 



-di - -a, - — ^{r) = EVL{r 



and 



if^$(0) = a2$(0) 

where and Qa,i3 are as in (12.61) and (12.71) respectively. 

The solutions to the radial equation can be obtained by taking 

n(r) = r^e-"'"'/2F(r). 

The resulting equation for F is 



-dl +\2ur + 



(27 + 1) 



p 2_ 2 

dr+ \ ^ +2u;(7 + l) 



Choose 



and set 



to find 



7 = k\a\ 

E = S + 2uj{k\a\ + l) 



{2k\a\ + 1) 



SF. 



Performing the change of variable 



y = ujr, 



(13. 9 p becomes 



yd^F + (1 + k\a\ - y)dyF + ^F = 



(3.2) 
(3.3) 
(3.4) 

(3.5) 

F = EF. (3.6) 

(3.7) 
(3.8) 
(3.9) 

(3.10) 
(3.11) 
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which is recognized [22j as the equation that has the Laguerre polynomials Lm"' as 
regular solutions if 

£ = Aum. (3.12) 

The radial part of the wave functions which are continuous at the origin for cu > and 
vanishing at infinity are thus given by 

n{r) = M^Y^^\y), (3.13) 

where 

Y::\'^\{y) = {yf\'^\''e-y/'Ly{y) (3.14) 

and Mm is a normalization constant that we shall specify later. From (13. 8 p and (I3.12p 
the eigenvalues are 

Em = 2uj[2m + k\a\ + I]. (3.15) 

Let us know turn to the angular equation (13. 4p . This equation has been solved in 
[19] by obtaining the eigenfunctions $ of Qa^p '■ 

g„,^<l>(0) =a$(0). (3.16) 

Clearly, H^<^ = a^$. One first observes that 

$0 = iVol sin 0|"/2cos'^/2 ^ gin 0)1/2 ^3_^7^) 

satisfies 

Qa,/3$o = -^(« + /3 + l)$o. (3.18) 
One then carries out the "gauge transformation" of Qa,i3 with $o: 

Qa,p = $0 'ga,/3$o; $ = %G (3.19) 
and after performing the change of variable 

x = sin0 (3.20) 
one finds that the functions G obey the equation 

Qc^^lsG = aG (3.21) 

with 

Q^,p = (1 - x)d,R - ^(1 - i?) - ^±1±1r. (3.22) 

This equation fl3.22p is recognized to be the differential-difference equation [151 [12] that 
has the little -1 Jacobi polynomials Pn"'^\x) as regular solutions. The eigenvalues are 

-n - ^^±f±i n even 

nodd • ^^-^^^ 

These polynomials have the following expressions in terms of hypergeometric 
terminating series 



[X) = Xn 



_n n+a+/3+2 \ / 1 _ ™ n+a+/3+2 



(3.24) 
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for n even; and 

{a + I3 + I)x 



l~n n+a+/3+l 
2 2 



a+1 i ^ 



a + l 



1-n n+a+/3+3 
2 2 . ™2 

2 



(3 



for n odd. 

Recall that the series 



2^1 



a, b 



n=0 



n!(c)„ 



(3.26) 



terminates whenever a or 6 is a negative integer. Here (a)„ = a{a + 1) ■ ■ ■ {a + n — 1) 
etc. are shifted factorials. 

The coefficients Xn are chosen so as to make the polynomials Pn"'^\x) monic, i.e. 



x"^ + 0{n — 1). They are given by 
(-1)5. 



Xn 



a + l ' 



.t + f + f + l)n 
/ a + l ' 



-1) 



' li + ^ + li- 

1+1 \ 

Ti+1 V 2 J H+i 

(^^ + f + 1 + l) n+1 



n even 



n odd 



(3.27) 



For a > — 1 and /3 > —1, the polynomial P^'^\x) are orthogonal with respect to the 
weight function 

uj{x) = |xr(l - x2)(''+i)/2(l + x). (3.28) 
The solutions to the angular equation (13. 4p are hence given by 



with 



TV 

$„(0) = -fX„(sin0), 

JVn 



X„(sin0) = $0(0;) Pi"''') (sin, 



(3.29) 



(3.30) 



and where Nn are also normalization constants that we shall specify below. Finally, the 
eigensolutions of (13. ip are 



^ r 



iVo 



(3.3i; 



iVo 

and the corresponding spectra are 



2uj 



2m + k \ n + 



sin0|"/2 cos^/2 0(l + Sin0)l/22^^a„|^^^2^)p^a,/3)(gi^ 

a + /3 + 1 



+ 1 



(3.32) 



Very much like the in the case of the TTW systems, the energy eigenvalues are linear 
in the quantum numbers n and m. For rational k, there is a degeneracy of states that 
we shall explain in the next section. 

The inner product on the space of wave functions is naturally defined by 



{fir,e),gir,e)) 



rdr 



d<P fir, <P)g{.r, 0). 



(3.33) 
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It can then be checked that the eigenf unctions of H^. (with k fixed) are orthonormal 



M^M^,Nn'Nn 



(3.34) 



where 



(X„(x),X„.(x)) 



X] 



(3.35) 



is the inner product on the space of wave functions of H^. Owing to the orthogonahty of 
the polynomials P"'^(a;) with respect to the weight function (13.281) . the functions Xjyx) 
are orthogonal. The norms have been evaluated in [19j. Recall x = sin*^, and we have 



iXn{x),Xn'ix)) = / d(sin0)|sin0|"(cos0)^"^(l + sin0)P°'^(sin0)P„";^(sin( 

/I 
dx Ixl^il - x^)^ {1 + x)p:^^{x)p^;^{x) 



5 ^ 

On,n' ^2 
n 



The normalization constants are given by [19j : 



with 



JVo(f + f + lJ 



7Vo(f + f+l) 



for n even, 
for n odd, 



2 2 /A 2 ' 2 



iVn 



r(f + f + i) 



1/2 



r(f + i))r(f + i)) 

and where T{x) is the standard gamma function. As a result 



n,m) 



N^Nn'MmMm' 



2u 



dyY!:\^-\{y)Y'S-'\y) (X„(x), X„,(x)) 
dyY!:\'^-\{y)Y'^r"\y) 
dy (l/)'l'^"le-^L^-l(2/)L^-l(2/) 



n,n' ^m,m' 



by the orthogonality of the Laguerre polynomials and assuming that 



2ujm\ 



(3.36) 



(3.37) 



r(m + /c|a„| + 1) ' 

where we have appended the index n to show the dependence on both m and n. 



(3.38) 

(3.39) 

(3.40) 

(3.41) 
(3.42) 

(3.43) 
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4. Superintegrability 

We now want to prove that for rational values of the Hamiltonians Hk are maximally 
superintegrable in addition to being exactly solvable. To that end, we shall explicitly 
identify the additional constants of motion that are at the root of the state degeneracies. 
We shall also determine the defining relations of the symmetry algebra formed by these 
conserved quantities. The supercharge Qa,i3 will be found to play an important role. To 
proceed, we shall use the recurrence relation approach that has been developed in [5]. 
This method uses the structural properties possessed by the separated eigenf unctions, 
to obtain the constants of motion from appropriate combinations of the corresponding 
ladder operators. For the Laguerre polynomials, these relations are known p3], as for 
the little -1 Jacobi polynomials, they have been determined in |18j . 

We look for recurrence relations on the wave functions that leave the energy fixed, 
i.e. automorphisms on the energy eigenspaces. To this end, assume k = p/q, then the 
transformations 



m ^ m — p, n ^ n + 2q 

do not change the energy: 
2u 



E 



{2mq + pn + p{a + /3 + 1) + g) 



(4.1) 
(4.2) 

(4.3) 



Recall, that the wave functions are indexed by m, n as 



m 







iy)Xnix) 



y = ojr , X 



with 



>^™'""'(l/) = (z/)"""'/'e-^/2L^""l(y) 
X„(x) = $o(x)Pi"'«(x). 



4.I. Ladder operators for the radial part 

Based on the ladder operators for the Laguerre polynomials 
ladder operators for the functions Ym"'"\y) 



i^_,|,„l,^y^l'^"l(l/) 



where E = Er, 



E 



1 



:i + kK\)dy-^^ 2y 



A;|a„|(l + k\an\) 



1 - k\an\)dy 



E 1 

4^ + - 



{m + l){m + k\an\)Y^ 



fc|a,i|— 2 
m+1 



we have the following 



r^i-i(y) 



rj;i-i(y) 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



= 2a; [2m + fc|a„| + 1]. To obtain the desired shift in the quantum 
numbers (14. ip we take a composition of K^s p times, with the corresponding value of 
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k\an\ shifted in each successive apphcation 

^l\a^\,E = Kk\a„\+2{p-l),E " " " -K'fc|a„|+2,£;-K'fc|a„|,£; (4.8) 
^-k\an\,E = K-(k\an\~2ip-l)),E " " " K ^(k\a„\-2),EK -k\a„\,E ■ (4.9) 

The action of these operators on the wave functions is given by 

= (-i)^i;^L";i+^^x„ (4.10) 

^-fc|a„|,E^r^'""'^n = (-'^Ti^ + l)p(^l«nl +^-P+ l)p V;t;'~''X„.(4. 1 1) 

It is important to note that ahhough the quantity is a function of m and n, it is 
unchanged by the operation m — )■ mil, and /c|an| — )• fc|an| =F 2. Hence, the change does 
not need to be accounted for in the successive apphcations of the operator. 

4-2. Ladder operators for the angular part 

Next, we look for operators which will change the value n by ±2g. To obtain these, 
we use the ladder operators of the little —1 Jacobi polynomials, as given in ^18j. The 
operators are: 

J+ = (x + (1 - x)R) (2Qa,i3 + l) + a + /3 

J_ = (x - (1 - x)R) (2Qa,(3 - l) + a - /3 (4.12) 

It is straightforward to verify that these operators satisfy 

{ J+, Q«,4 = - J+, { J_, Q^,^} = J_ (4.13) 

and so the operators J_|_ and J_ fix the space of P„ but change the degree of the 
polynomials by plus or minus 1, depending on the parity of the degree, n. The action 
of these operators on the polynomials P„ is 



2n(a+/3+n) 



Pn-1 for n even. 



^ 2(a + /3 + 2n + 2)P„+i for n odd, ^ ^ 

-2(a + P + 2n + 2)P„+i for n even, 

J^Pn - i 2(.+n)(/3+..) p ^ . (4.15) 
(a+p+2n) " ^ 



The operator which will change the degree of X„ by 2q is given by 

CnXn-2q for n CVCU, 
dnXn+2q for H odd,. 



where 



(4.16) 



(4.17) 



2g 



(i„ = (-2^)^ + n + 1^ . (4.18) 

Writing the J's in the the opposite order, we obtain 

J-'- ^ 4„(J,J.)'t-.Y„ = I "'"-f""' f " •'"2' (4.19) 

c„X„_2q for n odd. 



An infinite family of superintegrable Hamiltonians with reflection in the plane 



10 



where 

On = ' ' J 7^ ' ' '\ (4.20) 

dn = dn. (4.21) 

Thus, we have operators with the desired action on the basis states, however, their 
action depends on the parity of the quantum number n. 

4-3. Parameter independent ladder operators 

Next, we would like to combine the ladder operators for the radial and angular parts 
to obtain operators which give automorphism of the energy eigenspace and which are 
independent of the indices m and n. On the one hande, note that while the action of 
the operators 7=*='^^ depend on the parity n, the operators do not. On the other hand, 
the operators -f^'±fe|a„| e depend on the values of n and m. To remove this dependence, 
we shall replace |a„| by the operator Qa,i3- This has the added advantage of making the 
action of the K"^ depend on the parity of n. Thus, let us define new operators which 
are now differential operators in x and y with a reflection in x: 



KkQ,E 



r fc I QjTL 



K.^..YS-^X. = I f " (4.22) 

I Kk\ar,iEym 'a„, for u odd. 

The repeated application of this operator, denoted K^q ^, is defined as 

^kQ,E = KkQ+2(p-l),E ■ ■ ■ KkQ+2,EKkQ,E (4.23) 

and will act on the basis as 



TCP yfclani X - ! ^^-tknl,^^-'""'^" for n even, 

^kQ,Eym ^P^^^^^y;^K|^^ for U odd. 



(4.24) 



Similarly, 



T^v ^^k\a\^ f -^fela lE^^'""'^" for n even, 
I ^-fe|a„|,£^- for n odd. 

Note that the operator KkQ,E still has m and n dependence via the energy E and so the 
next step is to remove the energy from the operator. To do this, we push the constant 
E to the right and replace it with the Hamiltonian as 

KIq,h = KIq,e\e=h. K^_kQ,H = KIq^e\e=h. . (4.26) 
This operator, K^q ^ is parameter independent and acts on the basis by 

^^Q,if>;^'""'^n = KlQ^Ey^'^-\Xr, (4.27) 

^ r K\^^ |^F^''^"lx„ = (-l)P(m + l)p(A;|an|+m-p+l)pF^^;l"'^X, for n even, 
I = i-^n!:i-f"'Xn for n odd. 
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Hence, we have parameter independent operators Hi and H2 given by 

Si ^ %{J^J_Y%'KIq^^, ^ <!>o{J-J+r^o'K'-kQ,H, (4-28) 



which act on the basis as 

rfc|a„|-2p 



" ""'X^lX-t^'^n^-. for -Odd, ^'-''^ 

S.y^l-IX„ = ( " (4.30) 

I Cn ^n-2g for n odd. 



■ m+p 

with 



-(4.31) 



2q 



«,„ = 4,„ = (-in-Z")' ( + n + 1 ) . (4.33) 



"'2q 



2q 

From the action on the wave functions, we can use a standard Wronskian argument as 
in [5] to show that the operators Hi and H2 commute with the Hamiltonians. However, 
it is also possible to prove this explicitly as we have done in the appendix. 

Also from the action on the basis, we can verify that the operators Si and S2 are 
mutual adjoints. Since the action of Si and H2 on the basis involve only one term, it is 
enough to show that, for n even 

,n—2qi ra,n ) — {^2^ m+p,n-2q, ^n^,„) (4.34) 

and, for n odd, 

) — {^2^ m-p,ri+2q^ ). (4.35) 
For n even, the requirement fl4.34|) is equivalent to 



n m,n m+p,n—2q 



N"^ /\/f2 

n—2q m+p,n—2q 



(4.36) 



This can be directly verified since 

N'M^ {-{^)-n)l 

n m,n v \ 2 / / 2g 



iV^2.M^4-p,n-2, (-1) q (-^), (^), {^)^ (m + lUm - p + k\a^\ + 1), 

_ ^~^m+p,n-2q r^^-^ 



For n odd, the requirement fl4.35p is equivalent to 
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Which again is directly verified as 

(4.39) 



'2g 

"7+ ' 



Thus, Sj = S2 and we have a superintegrable system with if, _g and self-adjoint 
operators Si + S2 or — H2). 

y^.y^. Algebra Relations 

From the expansion coefficients (I4.31l) - fl4.33l) . we can use the action of the operators on 
the basis Ym''"^Xn to determine the algebra satisfied by the operators Si,S2 and Qa,/3- 
The commutator of Si with Qa,i3 is given by 



Similarly, 



-lanl^lYm' — Qa,l3^m,n^m+p ^ Xn-2q for U eVCU, 

I'^nl'^'l^m' "'^n — Q a,l3(^m,rXrn-'p^ ^ Xn+2q for U odd, 

(-|a„| + |a„_2g|)Cn^m+p ^^''-^n-2<? for u eveu, 

(knl - |an+2g|)Cn^m-p "^^^^^+29 for H odd, 

2gSi^^,„. (4.40) 



|a„|H2Fm'""'X„ - QaMm n^m-p ^^^-^n+2<? for U eveU, 



|o„|^l^rn''*"'^n ~ Qa.^S^^m+p' ^^^n-2i3 for U odd, 

(-|a„| + |a„+2g|)^+^„Fj;'!;'+^^X„+2g for n even. 



DCnWp ^^n-2g for U odd, 

= 2gS2^„^,n. (4.41) 
From the action on the basis, we can conclude that 

[Si, QaA = -2gSi, [S2, QaA = 2gS2. (4.42) 

Finally, writing the expansion coefficients in terms of e = m + so that E = 
2uj (2e + Mfij^ + 1 j , they become 

fcn ^ / + « + + \ (-2^)^(-f),(-^),(^),(^) 



^,n = {e-Y + ^)p [ 2 ^ 

, kn . f k{n + a + (3 + l) \ 
.„ = (e-^ + l)p ^ ^ ej 



2 "V 2 

fa,,, / k{„ + a + fl+l) \ (-2'')'(V),(i"¥^),(-^),(-^), 



P V 2 ■V2g 



-iy(-2y{^^ + n + l] . (4.43) 



2g 

Using these forms of the coefficients fl4.43p - fl4.43p . the action of [Hi, S2] on the basis is: 

iz: ]\tj — j i^m,n^m-p,n+2q ~ ^m+p,n~2q^m,n)^ m,n fo^ ^ eVCU, 

[^1, ^2j m,n - S „+ for n odd ^ ^ ^ 

L \^m,,n'^m+p,n-2q '-m-p,n+2q'-m,n) ^ rn,n iui li uuu. 



An infinite family of superintegrable Hamiltonians with reflection in the plane 13 
where 



"m,n m—p,n+2q 



x('l + !i') fl+° + '^ + "'l + f° + "-' I , (4.45) 



"m+pjn— 2g m,n 



/ q \ / q 

k{n + a + /3 + 1) 



2 



— e 



A-a-^\ (l-P-n. ^^^g^ 



2/a V 2 



q \ / g \ / 9 



'-If 2««(e - ^ + 1). (^tt±^±l±}l _ , 



^m+'p,n—2q^m,n V / V 2 



X ( i^)^ ( ^-";^-^ )^ (-^)^ (-^)^, (4.47) 
-lf2«Ml-e + -^ / + « + + 



p \ / p 



\ I -r ^ ' 1 + 1 + . (4.48) 



2/V 2 yv 2,, 

/g\ /g^ /g 

These coefficients are polynomial in n and e and so can be written as polynomials in 
the operators Qa,i3 and E as 

_E k{a + (5 + l) 1 _ a + + l 

^~4J 4 2' ^ 

which is compatible with the substitution E = H and |a„| = —Qa,i3 for even and 
I On I = (5a,/3 for n odd. The commutators of the operators Hi and S2 can thus be 
expressed as 



^m,n^m-p,n+2q ^m+v.n-2a^m.n ) I g-^fc(c»+/3+i)-2i^ 2Q^^^+a+;3+i for U CVeU, 

/ ^- 4^ -2 (4.49) 

-1 for n odd. 



-m,n^m+p,n—2q ^m—p,n+2q^m,n j \ ^_ H-ujk{a+fl+i)-2oj „_ '^Qa,i3 ° Z^' 



To guarantee that these two forms of fl4.49p are equal, it is enough to verify that 
under the transformation ra— t-— (n + a + /3 + l) the coefficients transform as 

m,n m—p,n+2q m,n m+p,n—2q'> 

/" f~ (A ^0) 

^m+p,n—2q^m,n ' ^Tn—p,n+2q'^m,n' y-c.^^j 

This requirement ensures that replacing |a„| with ±Qa./3 will give the appropriate 
expansion (14.441) for both even and odd basis and so the two forms in (14. 49 p . Thus, we 
have shown that the commutator and analogously the anti-commutator is a polynomial, 
of degree at least 2p + 4q in the operator Qa,i3 and 2p in the Hamiltonian. 
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^.5. Case k = 1 

Let us present explicitly the results for the case k = 1. The ladder operators are given 
by 



{{x + Vl - xm) {2Qa,p + 2) + a + (^(x - Vl - x^) {2Qa,is - 1) + a - /3 



H 1 



(4.51) 



^2 



(^{x - Vl - xm) {2Q^^fi - 1) + a - (^(x + Vl - x2i?) (2g,,^ + 1) + a + 
X (1 - fcQ„,;3)a, - ^ + ^A;Q„,^(1 - Q,,^) . (4.52) 

Finally, we can compute explicitly the structure relations for the symmetry algebra 
for k = \: 



-2^1, \^2-,Qc 



2S 



2, 



(4.53) 



■i,S2] =48Q^_^-8 
«2 + /32 _ 3 



+ 2 



2 + (/3 - a - 3)(/3 + a - 3)(/3 - a + 3)(/3 + a + 3) - 3 ) g,,;^ 



Ql,, - 2^Ql, + ^{Hf - Au') (4.54) 



and 



{Hi,S2} = -8Q^. + 2 



if2 + 2a;2(^2 ^ ^2) _ 58^2 



Qi. + 8a(3Q 



(4.55) 



(2^2 + 2/32 - 22)//! + a;2((a2 _ /32)2 _ 50(^2 ^ ^2) ^ ^93^^ 



8a/3(g2 - 12u;2 
3IJ2 



(«2 _ /32)2 _ ^0(«2 ^ ^2) ^ g 



(ifi2_4c^). (4.56) 



Note that these algebra relations are degree 7 and 8 in the momenta, respectively. 



5. An infinite family of superintegrable systems with 1/r 

As was shown for the TTW system [18], coupling constant metamorphosis 
be used to map this new infinite family of superintegrable systems Hk (12. ip to 
family of superintegrable systems Hk with a Coulomb term, 

Hu = -dl--d, + ^- + ^H^. 



p 



p 4p2 2 



To see this mapping, take the Schrodinger equation 

U2 



1 A;2 
-d^. - -dr + uh^ + -p^H^ 1 ^ - = 



and solve for 072/4 as 



E fc2 



B6] can 
a different 



(5.1) 



(5.2) 



(5.3) 
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Then, making the change of variables 

= p, (5.4) 
and interchanging the couphng constants with the energies 

E = -47, = E, (5.5) 

we obtain the Hamiltonian (15. ip with Schrodinger equation 

Hk^ = E^. (5.6) 

The wave functions in (15. 2p and (15. 6 p differ only in the change of variables (15. 4p and 
parameters (15. 5p . The separation constant A;|a„| becomes fc/2|a„| and so the quantization 
of the spectra E = u{2m + k\an\ + 1) becomes 

/ — - k 

-47 = 2V-E(2m + -|a„| + l) (5.7) 

or 

872 

(4m + A;|a„| + 2)2- 
The wave functions \E' are given by 



Furthermore, as was shown in [THl |26] the integrals of motion are preserved under 
this transformation. Alternatively, we can use the recurrence relation method described 
above to obtain the integrals, since the wave functions are essentially the same. Note 
that the integrals of iffc will be mapped to integrals of 11^12- 



6. Conclusion 

In this paper, we have constructed an infinite family of exactly solvable Hamiltonians 
indexed by a parameter k whose wave functions can be written in terms of Laguerre and 
little -1 Jacobi polynomials. For rational fc, we used the recurrence relation approach 
of Kalnins, Kress and Miller |5] to construct the integrals of motion out of the ladder 
operators of the two families of orthogonal polynomials. This method not only gives 
explicit expressions for the integrals of motion but also gives the symmetry algebra 
which they satisfy. 

We have thus provided one class of Hamiltonians with reflections corresponding to 
superintegrable systems in the plane. It would be of interest to look systematically 
for analogous models, that is to search for all separable Hamiltonians that are 
superintegrable when the occurrence of reflection operators is allowed. This would lead 
to the identification of orthogonal polynomials that fall outside of the Bochner Theorem 
in dZ]. 
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Appendix A. Direct proof of the commutation relations 

Let us prove explicitly that the operator Si commutes with if^. The proof for S2 is 
analogous. 

Theorem 1 The operator Si {4-28) commutes with the Hamiltonian Hk ^2. l\l for 



k = p/q. 

To this end, we prove by induction the following lemmas: 
Lemma 1 The following identity holds for all i : 



y 

Proof: This is a proof by induction. Recall the definition f l4.19p 

J+'2^ = %iJ+J^)%\ (A.l) 
As a direct result of the anti-commutation relations of the operators J± fl4.13p . we obtain 
[j+'',Q^,p]=2J+'' (A.2) 
[J^'\Ql,]=-^J'''\l-Qa,,). (A.3) 
Next, we use the form of the Hamiltonian (12. ip in terms of the variable y = ur"^ : 

Hu = uj [-Aydl -Ady + y+ ■ (A.4) 

And so, since J± have no dependence on y, their commutator with Hk is given by 

[J-'''',H,] = ^[J^'^\Ql,]. (A5) 
The final step is then to prove by induction that 

[j^^'',Ql,]=-U^^''{f-iQa,,). (A.6) 

The case £ = 1 is shown in (IA.3P and so we assume (IA.6P for i — 1 and compute 



= - 4J+'^^(l - Q„,,) - 4J+'2^-2 ((£ - If -{I- 1)Q„,^) J+'2 (A.8) 
= -U^^^\e-iQo.3)- (A.9) 
Thus, by induction we have proved the lemma. 

□ 
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Next, we will prove an identity for the commutator of the operator K^q ^ with Hj^. 
Along the way, we will give an explicit factorized form for the operator K^q ^. 



Lemma 2 The following identity holds for all I : 

Proof: First, we note that for £ = 1, the operator KIq ^ is 



^kQ,H 



H 1 



(A. 10) 



(A.ll) 



It is then straightforward to verify, using the explicit forms of the operators (IA.4I) and 
f lATTl) . that flAlOl) holds for £ = 1 : 



We shall prove the lemma again using induction, so assume ( lA.lOl) for 
purposes of induction, we take the factorized form of K^q^j as in f l4.22p 



^kQ,H — KkQ+2l^2,HKlQjj 



(A.12) 
1. For the 

(A.13) 



where KkQ+2e~2,H is defined as the operator with action on the basis as 

K,Q+2i^2,HKiQ]j,^ = K,Q+2i-2,EKiQ]E^, (A. 14) 

i.e. we move the constant E to the right and then make the replacement with H. 
Explicitly, we have 



kQ+2l-2,HK^Q^^ 



+ 



[^kQ,H^ Hk] 

Hk 1 



(A.15) 



{kQa,i3 + 2i- l)dy - ^ - ^ikQa,p + 2£- 2){kQ^^p + 2i-l) 



^^kQ,H- 



(A.16) 

If we use the induction assumption, (lA.lOp . then the final term \K^j^^,H}?\^ is 
proportional to K^q^^ and we obtain an explicit expression for KkQ+2i-2,H as 



Hk k^Qlfs - - 5)Qa,/3 - 2{3£ -m- 1) 



KkQ+2e~-2,H = {kQa,i3+2i-l)dy- 2^ 
Using this expression for KkQ+2£~2,H, we directly compute 

J^kQ+2l-2,H — J^kQ+2e~2,H 



(A.17) 



[KkQ+2e-2,H, Hk] = Au 
which in turn shows 



{i-l)kQa,^ 



y 



(A.18) 



[^kQ,Hi Hk] - KkQ+2e-2,H [kIqIh, Hk] + [KkQ+2e-2,H, Hk] K^^^^ 



Au{e + (ikQ) . 



kQ,H- 



(A.19) 



and so we have proven by induction that the expression (lA.lOp holds. 



□ 
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Proof of Theorem 1 The proof of theorem immediately follows from the two 
lemmas, li k = p/q then, 



4a;(p2 + fg) 



r^P _ j+,2q ( W(g'-gQa,/3) \ 
^kQ,H \ q^y J ^kQ,H 



= 0. (A.20) 

□ 
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